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ABSTRACT 

A reflectance -based method was used to provide an analysis of the in-flight 
radiometric performance of AVIRIS . Field spectral reflectance measurements of the 
surface and extinction measurements of the atmosphere using solar radiation were 
used as input to atmospheric radiative transfer calculations. Five separate codes 
were used in the analysis. Four include multiple scattering, and the computed radi- 
ances from these for flight conditions were in good agreement. Code -generated 
radiances were compared with AVIRIS-predicted radiances based on two laboratory 
calibrations (pre- and post-season of flight) for a uniform highly reflecting 
natural dry lake target. For one spectrometer (C) , the pre- and post-season cali- 
bration factors were found to give identical results, and to be in agreement with 
the atmospheric models that include multiple scattering. This positive result 
validates the field and laboratory calibration technique. Results for the other 
spectrometers (A, B and D) were widely at variance with the models no matter which 
calibration factors were used. Potential causes of these discrepancies are 
discussed. 


1 . INTRODUCTION 

We present here an updated analysis of the in-flight radiometric performance 
of the Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) 1 ’ 2 . AVIRIS is a 
scanning imager (400 - 2450 nm) that employs four spectrometers connected to common 
foreoptics by optical fibers. The spectrometers employ linear detector arrays and 
are labeled A (400-710 nm with 32 detectors), B (680-1280 nm, 64 detectors), C 
(1240-1860 nm, 64 detectors), and D (1830-2450 nm, 64 detectors). A test flight was 
carried out in mid-September 1987, over a field site at Rogers Dry Lake, California 
(Edwards Air Force Base) with the following objectives: (1) develop radiometric and 
spectral calibrations for the instrument; (2) determine in-flight signal/noise 
characteristics; (3) compare various methods of reducing the data to ground 
reflectance. To. support these studies, the following field measurements were made 
simultaneously with the overflight: (1) atmospheric optical depth using solar 
radiometers; (2) total precipitable water (in cm) with a spectral hygrometer; (3) 


SPIE Vol. 924 Recent Advances in Sensors. Radiometry, and Data Processing for Remote Sensing (1988} / 179 

197 


the ratio of direct to diffuse radiation incident at the surface using a hand-held 
ratioing radiometer (HHRR 3 ) . Spectral reflectance factor measurements were made of 
the bright playa surface using the Portable Instantaneous Display and Analysis 
Spectrometer (PIDAS 4 ), a Barnes Modular Multispectral 8-Channel Radiometer (MMR ' ) 
and an Exotech Model 100-AX radiometer . 


To meet the first objective we used a surface reflectance-based method and 
compared the solar radiance at the instrument reflected from the playa target as 
determined from a laboratory calibration of AVIRIS , with the radiance predicted 
from that target according to an atmospheric model. Previously we used the LOWTRAN 
6 code 8 for these comparisons, recognizing by so doing that we could only account 
for a single order of scattering in the predicted radiance. In the present work, we 
employed the LOWTRAN 7 code , which provides a description of multiple scattering 
(with Rayleigh and aerosol contributions) together with gaseous absorption, we also 
used multiple scattering codes developed by Herman and Browning , and Diner and 
Martonchik 1 . The latter two codes also account for multiple scattering and 
absorption effects, but rely on user inputs to quantify these at each wavelength. 
No spectral gaseous absorption databases are included within the codes themselves. 


In addition to incorporation of improved scattering codes in the analysis, 
we compared the results of pre- and post- flight season laboratory calibrations of 
the instrument using the Rogers Lake data. Apart from sound procedure, the 
relevance of this comparison stems from the post-1987 flight season discovery of 
the detachment of optical fibers connecting the foreoptics to spectrometers A and B. 
The point in time of this occurrence during flight operations during the summer of 
1987 is not known. But these comparisons, together with a reanalysis of the 
laboratory radiometric calibration procedure itself, has helped to sort out 
explanations for the substantially reduced radiometric performance of AVIRIS 
previously reported 7 . 


We also present comparisons of independent atmospheric optical depth and 
surface reflectance measurements by the Jet Propulsion Laboratory (JPL) , University 
of Arizona (UA) , and United States Department of Agriculture (VSDA) groups. The 
field data obtained agree closely and serve to reduce or eliminate the possibility 
that discrepancies between the AVIRIS calibrations were due to errors in the field 
measurements . 


2. TEST SITE AND FIELD MEASUREMENTS 
2.1 Site, flight and image data 

Rogers Dry Lake (Lat. 34° 55 N; Long. 117° 50 W) is a playa located 104 
km north of Los Angeles, California., at an elevation of 692 m above sea level. The 
surface consists of silty mudstone which, from X-ray diffraction analysis, contains 
quartz, plagioclase feldspar, kaolinite, (abundant) analcime, and calcite . The 
visible and infrared reflectance is however dominated by the presence of hydrated 
iron oxides that are npt detectible in the X-ray diffraction data. During the 
summer and fall months, the playa is dry and the surface is patterned with narrow 
dessication cracks. Intermittent streams have dissected the eastern half exposing 
a series of varicolored brown and tan strata beneath the surface layer. Such 
erosion has lead to broad patchiness in the exposed surface. Nevertheless, large 
areas of relatively uniform surface properties remain; these constitute the standard 
targets employed in our studies. 
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AVIRIS overflew the test site on September 14, 1987, around 10:20 AM Pacific 
Standard Time (18:20 hr UT) at an altitude of 20 km above terrain. AVIRIS has a 
one-milliradian instantaneous f ield-of -view and a 30-degree scan angle 1 , which 
determines a 20 m pixel size and 10.5 km swath width of the ground. The test flight 
took place under cloud- free, but hazy atmosphere conditions. The calculated 
visibility was 26 km. 

The image data contain periodic noise, and at least three spatially distinct 
patterns are present. The most prominent of these has a period of about 10.3 lines 
per cycle (0.097 line" ) normal to the banding which trends at an angle of about 18° 
to the scan direction. A second band set has a period 12.8 lines per cycle (0,78 
line" ) at an angle of about 65 degrees; the third periodic variation is manifested 
within the second set with a period of one line per cycle (1 line’ ) along the 
flight direction. A complete analysis of the noise characteristics of the data will 
be presented elsewhere. 


2.2 Spectral reflectance measurements 

Coincident with the time of aircraft overflight, we measured the spectral 
reflectance of the playa target with PIDAS and with the broadband MMR and Exotech 
filter radiometers. PIDAS covers the spectral range 425-2500 nm with sampling 
intervals of 0.88 nm between 450 and 900 nm, and 4.7 nm between 900 and 2500 nm. 
The MMR has filters matched to Thematic Mapper (TM) band widths with central 
wavelengths of 486, 571, 661, 838, 1677, 2223 nm, and an additional band centered at 
1254 nm. 


The field sampling strategy employed by the UA and USDA groups has been 
described in detail by Slater et al. A rectangle - 4 x 16 pixels in size (20 m per 
pixel) — was marked on the surface, its length being parallel to the flight line. 
Cross patterns of 16 measurements were made within each pixel, and averages and the 
scatter in each were computed individually and for the total 1024 measurements . The 
field reference was a BaSO^ panel. These reflectance factor measurements were 
traceable to the National Bureau of Standards diffuse reflectance standard by 
calibration with respect to a pressed halon surface . With PIDAS, four 
measurements were made at the center of each pixel and combined to secure an average 
from 256 measurements. In practice, about 50 measurements were lost, and the 
averages were computed from 200 spectra. The field reflectance standard was a 
portable pressed halon surface prepared to NBS specifications. 

Figure 1 shows that there is good agreement between the determinations of 
the spectral reflectance of the playa target measured by PIDAS and by the MMR. The 
PIDAS reflectance data depict mean values ±1 standard deviation of the observations. 
The scatter results from at least two sources, real variations in the surface 
reflectance, and a variation in the number of shadowed cracks present in the 
spectrometer f ield-of-view from place to place. The gap in reflectance between 1800 
and 1900 nm arises from saturation of the atmospheric water band over that interval. 
The jagged spectrum seen beyond 1700 nm results from noise (since corrected) 
introduced by the second of two amplifiers covering the infrared spectral region. 
The higher value at 2220 nm returned by the MMR is more nearly correct judging by 
the general concordance between these values and the hemispherical directional 
reflectance of a surface sample determined with laboratory ratioing spectro- 
radiometers . 
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Figure 1 Surface reflectance of playa calibration target on September 14, 1987. 

Solid line is spectral reflectance given by PIDAS ±1 standard devia- 
tion, representing an average of 200 spectra. Triangles are reflect- 
ance determined over TM bandpasses with MMR, representing averages of 
about 1500 measurements . 


2.3 Measurement of optical depth 

Measurements of incident solar illumi nation at ten wavelengths (370, 400, 
440, 520, 610, 670, 780, 870, 940 and 1010 nm) were made with two Reagan-type solar 
radiometers. These observations were carried out between approximately 7:00 AM and 
12:30 PM Pacific Standard Time. Optical depths were derived independently from 
slopes of conventional Langley plots of the solar IR radiance versus the secant of 
the solar zenith angle. The data for the two Reagan instruments are shown in Figure 
2 and compared there to optical depths obtained from the transmittance functions of 
the standard mid- latitude summer (mis) LOWTRAN 6 model. 

3. LABORATORY RADIOMETRIC CALIBRATION OF AVIRIS 
3.1 Definitions and procedure 

The laboratory radiometric calibration of AVIRIS was described in detail by 
Vane et al . The purpose of this calibrations was to obtain a set of functional 
relationships (one for each of 224 detectors) between radiance incident at the 
instrument, and instrumental output. If L is the radiance (in ^Wcm‘ nm' sr" ) and 
DR (digital response) the representation of instrument output, the relationship is 

L = DR/G (1) 

1 14 

where G represents a "gain" factor (G" = inverse responsivity ). In writing 
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Figure 2 A comparison of optical depths measured by JPL (•) and University of 
Arizona (A), September 14, 1987, Rogers Dry Lake. Optical depths at TM 
wavelengths (■) are calculated from the (A) data points. 

Eqn. (1) , it was assumed that detector dark responses were accounted for, i.e., zero 
instrument output for zero radiance input. 

The procedure for instrument calibration involved four steps: 1) radio- 

metric calibration of a spectroradiometer against a standard lamp, traceable to a 
National Bureau of Standards lamp; 2) production of light transfer curves (radiance 
vs. wavelength) for the AVIRIS radiance calibration source (40- inch diameter BaSO^- 
coated integrating sphere) using a calibrated spectroradiometer; 3) acquisition of 
AVIRIS DR output files as a function of integrating sphere radiance; and 4) 
generation of lookup tables of AVIRIS response vs. radiance for each detector at 
each wavelength. In practice, a linear relationship (Eqn. 1) was assumed and the 
lookup file generated from a single integrating sphere determination, thus 
determining a value of G. The radiometric calibration converted 10-bit DR values to 
32-bit floating point radiance values. 

We briefly recapitulate here the laboratory method used in (1) above because 
the equation employed by Vane et al . for calculation of spectral radiance seen by 
the standard spectroradiometer is incorrect. For its calibration the 
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spectroradiometer views a standard halon panel at an angle of 45 . The panel is 
illuminated at normal incidence by a standard irradiance source at a calibration 


distance above its surface. 


The irradiance E(CAL) at the standard panel is 

,2 


E ( CAL) - E(STD)(D std /D cal )‘ 


( 2 ) 


where E(STD) is the irradiance of the standard lamp provided by the manufacturer 
measured at the calibration distance Dg-jp. Assuming Lambertian diffuse reflectance 
from the halon panel with reflectance p, the radiance L(CAL) uniformly outward from 
the panel is 


L(CAL) - pE(CAL)/jr. (3) 

Combining Eqns. (2) and (3) gives the radiance of the spectroradiometer in terms of 
the irradiance of the standard lamp. This result is greater than that of Vane et 
al . 2 by a factor of 1/cos (45°), or Jl. The radiance of the integrating sphere Ljg 
is then determined in terms of the radiance emergent from the halon and the standard 
lamp irradiance from the relation, 

L IS " L(CAL)Ag R n SR ($g R (IS)/<&2 R (CAL))/A SR n SR ^ 

where $g R (IS) and $g^(CAL) are respectively the radiant fluxes observed by the 
spectroradiometer from the integrating sphere and from the calibration source. Ag R 
and flg R are the area of the entrance aperture of the spectroradiometer and its solid 
angular field-of-view. The etendue Ag R Og R is of course common to the spectro- 
radiometer for both the integrating sphere path and standard source path, but is 
included to emphasize that the determination is made from the total radiant flux 
received in each path. 

The absolute accuracy of the radiometric calibration was calculated to be 
about ±7% across the spectral interval 1800-2450 run, where the response of AVIRIS 
was observed to be most stable. Over the interval 400 to 1800 nm, additional 
uncertainties of a few percent are introduced because of thermal distortions 
(discussed below). 

3.2 Fre- and post-season instrument status and determination of G 

The initial laboratory calibration of AVIRIS was carried out in early June, 
1987, before the 1987 flight season. During this calibration, it was discovered 
that the output signal level from the four spectrometers varied with time in the 
presence of a stable input radiance. The instability was traced to thermally 
induced distortions from on-off cycling of heaters within the spectrometer barrels. 
The distortions produced a reduction in signal level without change in wavelength 
distribution of the energy. The magnitudes of these variations in signal output 
were: spectrometer A, 3.5%; spectrometer B, 7.7%; spectrometer C, 4.2%; 
spectrometer D, 2.4%. During subsequent engineering test flights, variations of 
twice these magnitudes were observed that also correlated with on-off cycling of the 
spectrometer heaters, and with distortion of the spectrometer base and instrument 
rack in-flight. 

After termination of the flight season (October, 1987), inspection revealed 
that the optical fibers connecting spectrometers A and B to the foreoptics had come 
loose from their mountings. Additionally, a longer wavelength blocking filter had 
fractured. The exact timing of the mishaps is not known. 
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A second post- flight season laboratory calibration of AVIRIS was carried out 
in October, 1987. In Figure 3, we compared the responsivities G (- DR/L) obtained 
from the pre- and post-season AVIRIS calibrations. The largest changes were seen to 
have occurred in spectrometers A and B with the post-season numbers smaller by 
factors varying between 0.36 and 0.91, compared to pre-season. Spectrometer C 
proved to be stable. The gain factor for spectrometer D decreased by a factor of 
about 0.8. 



Figure 3 Comparison of the pre- and post-season responsivities (gain factors G) 
for AVIRIS determined from the laboratory calibrations. 


4. COMPARISON BETWEEN AVIRIS- AND MODEL- GENERATED SOLAR RADIANCES 

An Important test of the validity of the laboratory radiometric calibration 
for flight conditions is the calculation of solar radiance at the instrument 
reflected from a target of known reflectance, when account has been taken of the 
atmosphere, as compared to the radiance predicted from the instrument response. In 
this section, we present: (1) the responses implied by the two AVIRIS laboratory 
calibrations; (2) the results of atmospheric modelling using the measured surface 
reflectances and atmospheric optical depths; and (3) a comparison of these various 
results . 


4.1 Radiances predicted from AVIRIS using pre- and post-season 
calibrations 

In Figure 4 are presented the spectral radiances produced for the playa 
target from the pre- and post-season laboratory radiometric calibrations. For each 
the ±1 standard deviation lines refer to scatter introduced because of natural 
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Figure 4 Spectral radiances predicted for the Rogers Day Lake target from AVIRIS 
from the pre- and post-season radiometric calibrations. 


variations in reflectance over the 64 pixels used in the sample and to instrumental 
noise. As implied by the changes in responsivity shown in Figure 3, the post-season 
calculated radiance for this target was greater over the spectral regions covered by 
spectrometers A, B and D, and essentially unchanged in the region of C. 

4.2 Radiances predicted by atmospheric models 

We have employed four atmospheric radiative transfer (RT) codes in these 
studies. The results from each are compared below: 

g 

(1) LOWTRAN 6 calculates atmospheric transmittance and radiance for a given 
atmospheric path. The code includes molecular absorption and single scattering from 
aerosols and molecules. The spectral resolution of 20 cm’ is sufficiently greater 
than the spectral sampling interval of AVIRIS (10 nm) to provide a very useful basis 
for wavelength calibration and effective spectral response for the instrument under 
flight conditions. The standard mid-latitude summer model was employed, and the 
geographic and solar conditions were the same as for other codes used, as given in 
Table 1. 

Single scattering was accounted for using the following equation, 

^total ~ ^"scattered + ^reflected 

where L scatterec j is the radiance from single scattering in the atmosphere above any 
ground elevation, and L re fp ecte d is the ground- ref lected direct radiance reaching 
the sensor. Two problems exist with the single scattering option in LOWTRAN 6. 


186 / SPIE Vol. 924 Recent Advances in Sensors, Radiometry, and Data Processing for Remote Sensing (1 988) 

204 



First, only a single target reflectance is accepted at all wavelengths in the 
unmodified code. We modified the program to allow for ground reflectance that 
varied with wavelength as given by PIDAS. Second, the program does not compute 
ground reflected radiance unless the target is at sea level. To circumvent this 
difficulty, the downward and upward transmitted parts from the target elevation were 
computed separately, and multiplied by the surface reflectance. The equation is 

^reflected ” ^O^o/^^dn^up^PIDAS 

where E Q — exoatmospheric solar irradiance, — cosine of the solar zenith angle, 
T^ n — transmittance of downward path, T -* transmittance of upward path, and Fpi^S 
is the reflectance of the (horizontal) surface returned by PIDAS. 

(2) The principal improvement provided by the LOWTRAN 7 code, from our 
standpoint, is the inclusion of multiple scattering, although it has not yet proved 
possible to modify the program to account for ground reflectance at any elevation 
except sea level. This code has been changed to account for variable surface 
reflectance as given by PIDAS. 

(3) The UA RT code by Herman and Browning 10 has been used extensively in the in- 
flight calibration of Landsat providing internally consistent results to within 
2.5% . The code uses Mie theory to compute aerosol scattering functions for a 
specified size distribution, and the Gauss -Seidel iterative technique to solve the 
equation of radiative transfer. The aerosol phase function is sampled at discrete 
angular intervals. 

(4) The 5S code of Tanre et al. 16 was modified at the University of Arizona to 
allow user specification of terrain elevation and sensor altitude, including other 
changes. The 5S code makes extensive use of analytical expressions and preselected 
atmospheric models, resulting in a very short execution time. The code has proven 
to be significantly more accurate than models with fast or moderate execution times 
(Royer, et al . 17 ) . 

(5) The JPL code 11 also computes aerosol scattering functions from Mie theory 
and employs the Gauss-Seidel method, but traces spatial Fourier components of the 
radiation field, and represents the aerosol phase function as a 24 -term Legendre 
series. The JPL code is believed, in theory, to have the same degree of certainty as 
the UA code, but has only recently been tested against field data. 

The results of calculations for these five models are given in Figure 5. 
Vertical bars on the JPL model points represent approximately the variation in 
calculated radiance from the variation in surface reflectance as determined with 
PIDAS (Figure 1). The spectral radiance generated by LOWTRAN 6 and 7 were for the 
standard mid-latitude summer model . The spectral radiance generated by the JPL, UA 
and 5S codes were for the geographic, solar and atmospheric conditions summarized in 
Table 1. Note that the LOWTRAN 7 result refers to an atmospheric path from Sun-to- 
sea level to an altitude of 20 km. Removing the lower 690 m of atmosphere, 
corresponding to the actual elevation of Rogers Lake, can be expected to increase 
the computed radiance somewhat, because of decreased attenuation, but to decrease 
the scattered component because of the reduction in scattering volume. The net 
result might be expected to resemble more closely the Herman- Browning radiance, 
since this code accounts for molecular absorption through the 5S code, at the 
proper elevation. The LOWTRAN 6 result, which represents the sum of directly 
transmitted radiation reflected by the ground and singly scattered radiation in the 
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Figure 5 Comparison of radiances generated by atmospheric RT models: (1) 

LOWTRAN 6, single scattering; (2) LOWTRAN 7, multiple scattering; (3) 
UA code; (4) 5S code; and (5) JPL code. 


TABLE 1 


MODEL CONDITIONS EMPLOYED WITH UA, 5S AND JPL RT CODES 
FOR ROGERS DRY LAKE, CALIFORNIA SEPTEMBER 14. 1987 


SOLAR ZENITH ANGLE: 37.3 
SOLAR AZIMUTH ANGLE: 142.9 
SOLAR DISTANCE, AU: 1.0058 
LATITUDE: 34® 59.7' 
LONGITUDE: 1170 52.4' 
ELEVATION: 694 m 2276’ 
AIRCRAFT ALT.: 20 km 65620' 
TIME OF OVERPASS: 18:21 UT 
CALCULATED VISIBILITY: 26 km 


JUNGE SIZE DISTRIBUTION: 3.46 
AEROSOL SIZE RANGE: 0.02- 5.02 
REFRACTIVE INDEX: 1.54-0.011 
PRESSURE: 920.4 mb 
TEMPERATURE: 21»C 
H z O CONTENT: 0.729 gW/km 
INTEGRATED PATH H z 0: 1.76 g/cmZ 
OZONE CONTENT: 0376cm-atm 
NADIR VIEWING ANGLE: <5« 


THEMATIC MAPPER BAND 

1 

2 

3 

4 



5 

7 

CENTRAL WAVELENGTH (nm) 

486.3 

570.6 

660.7 

p™ 9 

1050* 

1250* 

1677.0 

2223.0 

OPTICAL DEPTH 

AEROSOL 

0.299 

0.237 

0.191 

0.135 

0.0958 

0.0706 

0.049 

0.032 

RAYLEIGH 

0.148 

0.077 

0.042 

0.016 

0.0065 

0.0032 

0.0010 

0.0003 

GAS TRANSMITTANCE, (5S) b 

0.988 

0.920 

0.962 

0.964 

0.999 

0.9965 

0.990 

0.979 

SURFACE REFLECTANCE 

0.253 

0.335 

0395 

0.430 

0.416 

0.426 

0.446 

0.409 

SOLAR IRRAO..Wnr 2 A ,nr 1 

1883.84 

1799.95 

1533.35 

1011.00 

658.19 

464.49 

221.13 

71.26 

NORMALIZED CODE RAD. 11 

.06797 

.07985 

.09622 

.1076 

0.1046 

0.1082 

.1120 

.1029 

RADIANCE, Wm-^m' 1 sf 1 

UA 1u 

126 3 

135.6 

140.6 

101.6 

66.9 

48.6 

24.2 

7.1 

5S 16 

128.7 

139.3 

144.1 

103.0 

67.2 

48.7 

24.3 

7.2 

JPL 11 

128.0 

143.7 

147.5 

1003 

68.69 

49.86 

24.2 

7.0 

PERCENT DIFFERENCE 

100 x (JPL - UAJ/JPL 

1.4 

5.6 

4.7 

-1.4 

2.6 

2.5 

0.0 

-1.4 

100 * (JPL - 5S)/JPL 

-.52 

3.8 

2.3 

-2.8 

2.2 

2.3 

-0.3 

-2.9 

100 x (5S - UA)/UA 

1.9 

2.7 

2.5 

1.4 

0.4 

0.2 

0.4 

1.4 


a MONOCHROMATIC VALUES 

b CALCULATED FROM A RT CODE OFTAURE etal. 16 USED HERE FOR DETERMINING GASEOUS TRANSMITTANCES 
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atmosphere is about 35-40% lower at shorter wavelengths than the results from the 
other codes. The reasons are twofold, LOUTRAN 6 ignores the diffuse irradiance of 
the sky on the target and the upward atmospheric path radiance. The distinction 
between the results for LOWTRAN 6 and LOUTRAN 7 essentially disappears at 
wavelengths greater than about 1300 nm. 

The JPL, UA, and 5S models have been exercised at the TM band wavelengths 
only. The numerical differences in output for these three models are summarized in 
Table 1. In general the differences are greatest in the region of strongest 
scattering and may arise in the method of describing the phase function. In par- 
ticular, for the Herman- Browning and JPL codes, the differences disappear for the 
case of a purely Rayleigh atmosphere ( J . Martonchik, personal comm.). Close agree- 
ment is also found between the continuum LOWTRAN 7 radiance and the model-computed 
radiance in the TM bands. A better overall fit between these envelopes could be 
achieved with with closer spacing of the computed points. Two additional points at 
1050 and 1250 nm have been added to realize this in part. 

4.3 Comparison of AVIRIS-generated and model-generated radiances 

The AVIRIS-generated radiance for the playa target using the pre- and post- 
season gain factors is compared to the three model-generated spectral radiance 
distributions in Figure 6. To be definite in the comparison, we will regard the UA, 
5S, and JPL codes as more accurately describing the radiance to be expected in TM 
bands 1-4 (spectrometers A and B) from a perfectly functioning and perfectly 
calibrated instrument because, as we have mentioned, the LOWTRAN 7 simulation does 
not account properly for the actual elevation of the site. In bands 5 and 7 
(spectrometers C and D) , atmospheric scattering is diminished as shown in Figure 5, 
and we can regard comparison with the LOWTRAN 7 model in a better light. The 



Figure 6 Comparison of RT code results with radiance for the playa target as 

determined from AVIRIS according to the pre- and post-season radio- 
metric calibrations. 
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differences between observed and the JPL- and UA-modelled radiances at TM band 
wavelengths, where these code radiances were unambiguously determined, are given in 


Table 2. 


TABLE 2 


Comparison (in %) of observed and modelled radiances at wavelengths 
of TM bands. (+) or (-) values designate observed radiance greater or 
less than % of the model radiance. (A) indicates TM band falls within 

region of spectrometer A, etc. 



Pre- 

-season 

Post- 

season 

Band 

JPL 

UA 

JPL 

UA 

1 (A) 

-43 

-40 

+47 

+52 

2 (A) 

-27 

-23 

+58 

+67 

3 (A) 

-32 

-27 

+40 

+47 

4 (B) 

-37 

-41 

+13 

+12 

5 (C) 

--15 

-- 15 

-- 15 

— 15 

7 (D) 

0 

0 

-+30 

-+30 


The comparison with LOWTRAN 7 is given in Figure 7. Here, for spectrometer 
C, the agreement between AVIRIS -derived and model radiance in the continuum region 
is within (-)5%, for both pre- and post-season calibrations. For spectrometer D, 
the pre-season result is about 25% higher and the post-season 50-75% higher than the 
model. The result for spectrometer C is encouraging and leads us to believe in the 
general correctness of the laboratory and model comparisons attempted here, despite 
unfavorable comparisons in the other spectrometers. 



400 600 800 1000 1200 1400 1600 1800 2000 2200 2400 2600 


WAVELENGTH, nm 

Figure 7 Percentage difference 100 (LOWTRAN -AVIRIS) /LOWTRAN for the playa target 

and pre- and post-season calibrations. Code used is LOWTRAN 7. 
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The region 450-1050 nm, covering all of the A and part of the B region is 
presented scale -expanded in Figure 8 and compared to the LOWTRAN 7 result which 
provides the spectral continuity required to interpret (in detail) the AVIRIS 
spectrum. Possibly excepting the concave upward portion of the observed AVIRIS 
spectrum between 450-550 nm, which may represent the Fe +3 band seen in the playa 
reflectance (Figure 1) , the additional features present were introduced by the 
instrument or are noise. Compare this result with that reported by Vane 18 (his 
Figure 6) at Stonewall Playa, Nevada, in which similar features are present. Note 
also the prominent discontinuity near 700 nm introduced between the spectral regions 
of A and B by the post-season calibration. 



Figure 8 Comparison of LOWTRAN 7 and AVIRIS -predicted radiances for the spectral 
region 450-1050 nm. The degraded LOWTRAN spectrum calculated from Eqn 
(7) using FWHM AA - 9.58 nm in Eqn (8). Post-season calibration. 


5. EFFECTIVE IN-FLIGHT SPECTRAL RESOLUTION 

We can use the AVIRIS spectra of Figures 6 and 8 to provide estimates of the 
effective in-flight spectral "resolution", and by inference the spectral sampling 
interval. By spectral resolution, we mean the detection of separate radiance minima 
associated with neighboring absorption bands, arising, in this case, in the 
atmosphere. The spectral sampling interval is obtained qualitatively by dividing 
the wavelength interval covered by each spectrometer by the number of detectors 
covering the interval. 
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An improved evaluation of the spectral sampling properties of AVIRIS can be 
obtained by examining the spectral response function. For each detector located 
provisionally at wavelength A^, the signal s(A^) generated by the incident spectrum 
L( A) is 

s(X t ) = f f(A £ - A)L(A)dA (7) 

o 

where f(A^ - A) is the spectral response or spread function of the system at A^ from 
the radiance at wavelength A. In practice, this is evaluated by restricting L(A) to 
a narrow band of wavelengths AA at successive wavelengths A to secure approximately 

s(A. - A') = f(X L - A')L(A')AA (8) 

The spectral response curves for a number of AVIRIS detectors have been determined 
by Vane et al. using AA ** 1 nm. The function 

f(x) - exp(-h 2 x 2 ) (9) 

proves to be a good representation for the normalized response s(A- — A)/L(A)AA 
(Figure 9) . The central wavelength of the channel is taken to be tne maximum of 
f(x), thus fixing A^, and the full width AA^ at half maximum height (FWHM) as the 
effective channel spectra width of sampling interval. This gives h — 1.665/AA^. 
The spectral sampling intervals for AVIRIS for all spectrometers have been 
determined by Vane et al. as follows: spectrometer A, 10.0 nm, spectrometer B, 9.58 
nm, spectrometer C, 9.86 nm; spectrometer D, 9.85 nm. 
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Figure 9 Comparison of laboratory-measured spectral response with approximation 

of the single-channel normalized spectral response by the function 
f(x) = exp[-h x 2 ] (smooth curve). 
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In principle, the resolution is determined from the sampling interval and 
Nyquist's sampling theorem. If the frequency of sampling is l/SA s nm" , then the 
sampling interval A s = 1/25A S . Thus, for 5A g = 10 nm, A g = 20 nm; that is, narrow 
bands in adjacent channels will produce a single minimum in radiance and are not 
distinguished unless separated by >20 nm. 

We now use these results to examine the implied resolution throughout the 
spectrum. From Figure 6, at 2050 nm (spectrometer D) , the pair of CO 2 bands are 
separated by 50 nm, and are clearly resolved, thus establishing a resolution <50 nm 
for that spectrometer. At 1600 nm (spectrometer C) , the atmospheric CO 2 bands there 
are separated by 30 nm and are also resolved. 

An improved description of the in-flight spectral resolution is possible 
taking account of the functional form of the spectral response curves, Eqn (9). For 
AVIRIS, a single detector sees a signal that is represented by Eqn (7), which can be 
used with Eqn (9) to calculate a degraded full-resolution LOWTRAN 7 equal to that of 
AVIRIS. By comparing the degraded LOWTRAN spectrum to the observed AVIRIS spectrum 
and using h as a variable parameter the effective AVIRIS spectral resolution can be 
deduced. This comparison is shown in Figure 8 using FWHM — 10 nm (h — 0.1665), 
which represents the expected AVIRIS spectrum if the instrument were actually 
sampling at the laboratory- determined estimated sampling interval. A better match 
(Figure 10) is provided by doubling the width of the Figure 9 spectral response 
curve so that the full sampling width at half -maximum is 20 nm. This is in accord 



Figure 10 Comparison of the observed AVIRIS spectrum and LOWTRAN 7 spectrum 
degraded using FWHM AA - 20 nm in Eqn (8). Post-season calibration. 
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with the post-season spectral calibration of AVIRIS, which revealed that, because of 
the detachment of optical fibers in spectrometers A and B and the resulting de- 
focussing, the spectral FWHM of these two spectrometers was 18-20 nm compared to 9- 
10 nm in the pre-season calibration. Spectrometers C and D were found to have 
remained unchanged. 


6. SUMMARY 

Data from an AVIRIS overflight of a field test site at Rogers Dry Lake, 
California have been analyzed to characterize the in-flight radiometric and spectral 
performance of the instrument. The combined field data sets of the UA, USDA and JPL 
groups, consisting of spectral reflectance observations of the playa and optical 
depth measurements of the atmosphere, are presented. Results of the measurement 
programs agreed closely (a few percent for reflectance, 1% for optical depth 
determinations). Results of two laboratory radiometric calibrations of AVIRIS, pre- 
and post-season are presented. Significant differences in gain factors were found 
between the two calibrations, varying from reductions of as much as a factor of 2.8 
in performance (spectrometer A) to essential stability (spectrometer C) . Employing 
a reflectance -based method, we compared the AVIRIS-generated response to five 
separate atmospheric RT models . Reasonable and encouraging agreement was found 
between the instrument-produced radiance using the post-season calibration data, and 
the three multiple scattering RT models employed, particularly for spectrometer C. 
A scenario that may explain much of the variable behavior of the other spectrometer 
is as follows: (1) instrument was calibrated pre-season; (2) optical fibers 
detached during flight season prior to September 14 calibration flight; (3) 
instrument was re-calibrated post-season with fibers detached but not necessarily in 
positions occupied during the September 14 calibration flight. While the fiber 
detachment problem may partly account for the variable performance of AVIRIS 
relative to that expected, a complete story must include the effects of in-flight 
thermal and mechanical distortion, since the flight environment is drastically 
different from that of the laboratory. There seems to be no way of separating these 
effects for the past data sets because AVIRIS has since been refurbished. We must 
therefore rely on future test flights to explore such questions. 
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